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Abstract 

Carnitine palmitoyltransferase IB (CPT1B) and adrenoceptor beta-3 (ADRB3) are critical regulators of fat metabolism. CPT1B 
transports free acyl groups into mitochondria for oxidation, and ADRB3 triggers lipolysis in adipocytes, and their respective 
polymorphisms E531K and W64R have been identified as indicators of obesity in population studies. It is therefore important 
to understand the effects of these mutations on ADRB3 and CPTIB function in adipose and skeletal muscle tissue, 
respectively. This study aimed to analyze the rate of lipolysis of plasma indicators (glycerol, free fatty acids, and beta 
hydroxybutyrate) and fat oxidation (through the non-protein respiratory quotient). These parameters were measured in 37 
participants during 30 min of aerobic exercise at approximately 62% of maximal oxygen uptake, followed by 30 min of 
recovery. During recovery, mean respiratory quotient values were higher in K allele carriers than in non-carriers, indicating 
low post-exercise fatty acid oxidation rates. No significant differences in lipolysis or lipid oxidation were observed between R 
and W allele carriers of ADRB3 at any time during the aerobic load. The substitution of glutamic acid at position 531 by 
lysine in the CPTIB protein decreases the mitochondrial beta-oxidation pathway, which increases the non-protein 
respiratory quotient value during recovery from exercise. This may contribute to weight gain or reduced weight-loss 
following exercise. 
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introduction 

A sedentary lifestyle has contributed to a substantial increase in 
the prevalence of overweight and obesity [1,2,3]. Although aerobic 
exercise can reduce excess body fat [4,5], genetic factors may 
determine the susceptibility of an individual to gain weight [6,7] 
and alter the beneficial effects of physical activity [8]. Aerobic 
exercise increases adipose tissue lipolysis, which releases free fatty 
acids (FFA) into the bloodstream, FFA uptake by active skeletal 
muscle, and FFA oxidation in the mitochondrial matrix [9]. 

Fatty acid oxidation during physical activity is strongly 
dependent on both the rate of lipolysis in the adipocytes [9, 1 0] 
and membrane translocation of long chain fatty acids (LCFA) into 
the mitochondrial matrix of the muscle cell [9,11]. Corpeleijn et 
al. demonstrated that intramyocellular lipid oxidation depends on 
extracellular FFA availability as well as mitochondrial function 
[12]. It has been documented, that FFA concentrations increase 
during the recovery period after exercise at 65% of maximal 
oxygen uptake (V02max), which releases more FFA for oxidation 



in skeletal muscle. This is supported by the decline in non-protein 
respiratory quotient (RQnp) value during recovery [13]. 

The beta-3 adrenoceptor {ADRB3) gene is expressed in visceral 
and subcutaneous adipose tissue and encodes the catecholamine 
receptor that activates lipolysis [14]. A polymorphic variant of 
ADRB3 is located on chromosome 8pl2 [15] with the substitution 
of tryptophan (W) by arginine (R) at position 64 {W64R). The R 
allele has been associated with excessive abdominal fat accumu- 
lation, higher body mass index (BMI), lower high density 
lipoprotein (HDL) levels, and insulin resistance [14,16,17]. 
However, in two population studies conducted on the effect of 
physical activity, the results were not different between W and R 
allele carriers [18,19]. 

Carnitine palmitoyltransferase 1 [CFTl; chromosome 
22ql3.33), expressed in skeletal muscle, heart, and the testis 
[20], regulates long-chain acyl group entry into the mitochondria 
for incorporation into the beta-oxidation pathway [9,20]. 
RobitaiUe et al. found a significant association between the 
heterozygous E531K allelic variant genotype, which replaces 
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glutamic acid (E) for lysine (K) at position 531 of muscle-type 
CPTl [CPTIB], and indicators of obesity when the proportion of 
dietary fat was S34.4%, indicating a gene-diet interaction [11]. 
Another study conducted with people of African desc:ent in 
Trinidad and Tobago found that K allele homozygote carriers had 
a lower percentage of intermuscular fat, but a higher percentage of 
subcutaneous fat when compared with heterozygous and homo- 
zygous E allele carriers [21]. Moreover, an association between the 
K allele and metabolic syndrome traits was found in an adult 
German population [22]. 

We hypothesized that under aerobic effort, R allele carriers of 
ADRB3 would present lower plasma concentrations of glycerol 
(Gly) and FFA, whereas K allele carriers of CPTIB would present 
higher values of RQnp. Indices of lipolysis and fat oxidation can 
also help determine whether the R and K alleles interact to reduce 
the fat oxidative capacity. This study aimed to determine the 
influence of these genetic variations R allele on lipolysis and that of 
the K allele on fatty acid oxidation by comparing carriers and non- 
carriers during and after aerobic load. Results indicate that the 
E5'ilK affect fatty acid oxidation regardless of W64R. 

Methods 

The Ethics and Biosafety Committee of Colima University 
reviewed and approved the study protocol in accordance with the 
Declaration of Helsinki. AH participants provided written informed 
consent. 

Participants 

In total, 45 participants were included from a group of 285 
individuals who had recently participated in a study of the 
association between markers of obesity and genetic polymor- 
phisms. Participants had normal hemoglobin, blood glucose, blood 
hpid, and hpoprotein levels, and fulfilled the additional inclusion 
criteria of BMI values and proportions of dietary macronutrients 
(carbohydrates, fats, and proteins) within the average range 
[mean ± 1 standard deviation (SD)] established in the association 
study (unpublished data). The proportions of dietary macronutri- 
ents were obtained by 24 h recall. Data were analyzed using Diet 
Analyst software (Parsons Technology Inc NUTRIDATA Soft- 
ware Corporation). We also calculated the quantities of saturated, 
monounsaturated, and polyunsaturated fatty acids (SFA, MUFA, 
and PUFA, respectively) in the diet adjusted to the body weight. 

Subcutaneous whole body fat was determined from sum of 9 
skin folds (^9P); this included triceps, subscapular, biceps, iliac 
crest, supraspinal, abdominal, front thigh, and medial calf in 
accordance with the guidelines and equipment of the International 
Society for the Advancement of Kinanthropometry (ISAK), [23] 
as well as forearm skinfold measured on the distal lateral line at 
maximum forearm girth [24]. Subcutaneous trunk fat was 
determined from the sum of 4 trunk skinfolds (24S); this included 
subscapular, iliac crest, supraspinal, and abdominal [23], and 
abdominal adiposity was determined by waist circumference (W C) 
measured in accordance with ISAK guidelines and equipment 
[23]. None of the participants had undertaken an exercise 
program during the past year, and none had received drug or 
hormonal treatments during the past 6 months. 

Evaluation of Cardiovascular Function 

Evaluation was performed by a cardiologist. All participants 
were apparently healthy; however, clinical examination revealed 
respiratory or cardiovascular anomalies in some participants, and 
these underwent further review with echocardiography and 



electrocardiography resulting in the exclusion of 3 participants 
who required treatment. 

Cardiopulmonary Fitness 

The Tecumseh step test was used to assess fitness on a 20.3-cm 
high step. According to its guidelines, participants performed 24 
steps/ min for 3 min and maintained a constant stepping rate with 
the aid of a metronome and assistant. Heart rate (HR) monitors 
(Polar-mod A-5 902805, Kempele, Finland) [25] recorded 
participant HRs prior to exercise, and in the recovery phase at 
1 and 3 min post-exercise. The HR value at 1 min was classified in 
accordance with established categories [26]. 

Genomic DNA 

The sodium dodecyl sulfate (SDS)-proteinase K method was 
used to extract genomic DNA from white blood cells, using a total 
of 7 mL of venous blood collected in two 5 mL tubes with EDTA 
(Ethylenediaminetetraacetic acid, 0.25 M). Briefly, cells were 
washed thrice in ammonium chloride (131 mM) plus ammonium 
bicarbonate (0.9 mM) with intervening centrifugation (15 min, 
5000 revolutions per minute (rpm), 4°C), digested in proteinase 
buffer (50 mM Tris, pH 7.5; 100 mM EDTA, pH 8.0; 400 mM 
sodium chloride (NaCl); 0.5% SDS; and 50 M-L of 20 mg/mL 
proteinase K), and agitated for 12 h at 37°C. To separate the 
DNA from protein, 6 M NaCl was added under centrifugation 
(15 min, 10000 rpm, 4°C). The supernatant (3 — 4 mL) was mixed 
with 15 mL absolute ethanol to condense the DNA. The DNA 
was transferred to 1 .5-mL Eppendorf tubes containing 1 mL of 
70% ethanol and then centrifuged for 10 min at 7500 rpm. 
Ethanol was removed by evaporation using a vacuum system, and 
the DNA pellet was resuspended in TE solution (10 mM Tris, 
pH 7.5; 1 mM EDTA) at 4°C for 12 h. DNA level was quantified 
by absorbance spectrophotometry at 260 nm. Aliquots (100 |i,L) at 
100 ng/|i,L DNA were stored at 4°C until used for genotyping. 

Genotyping of Polymorphisms 

The W64R and E531K were genotyped by polymerase chain 
reaction— restriction fragment length polymorphism (PGR- RFLP) 
analysis. The reference single nucleotide polymorphisms (refSNPs) 
for W64R and E531K are. rs4994 and rs470117, respectively 
[27,28]. Primers were designed using Primer 4.0 (http://frodo.wi. 
mit.edu) as foUows: F-AGCCAGGTGATTTGGGAGAC and R- 
GTCACACAGAGCACGTCCAC for amplification of a 429-base 
pair (bp) fragment including W64R and F-GCACTTGCCTCT- 
TACCCACA, and R-CTTCTTGATGAGGCCTTTGC for 
amplification of a 31 1-bp fragment including E531K, which were 
synthesized by Invitrogen Custom DNA Oligos (Life Technolo- 
gies). 

The hybridization temperature was 60°C for W64R and 63°C 
for E531K. Amplification was performed in 35 cycles, each 
comprising 30 s denaturation at 94°C, 30 s annealing at 60°C or 
63°C, and 30 s elongation at 72°C. The PGR product for W64R 
was digested for 16 h in 5 U PflMI enzyme (New England Bio 
Labs. R0509S. Lot: 0280905), which recognizes the tryptophan 
base of the W allele and cuts the product into two fragments of 252 
and 177 bp. After digestion, the full 429-bp fragment is therefore 
from the R allele. The PGR product for E531K v/as digested for 
16 h in 5 U Taql (New England Bio Labs R0149S. Lot: 0540804), 
which recognizes the glutamic acid base of the E allele and cuts the 
product into two fragments of 217 and 94 bp. Thus, the full 311- 
bp fragment is from the K allele. Digestion products were 
separated by 6% polyacrylamide gel electrophoresis and stained 
with 0.2% silver nitrate. 
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Blood Chemistry Analyses 

Venous blood samples were taken from the antecubital vein 
12 h after fasting [29], approximately 5-6 weeks before the 
aerobic load. The parameters measured were (mg/dL): total 
cholesterol (T-Chol), triacylglycerol (TAG), high density lipopro- 
tein cholesterol (HDL-C), and low density lipoproteins cholesterol 
(LDL-C) using standardized colorimetric methods (Spinreact SAU 
Ctra Santa Coloma, 7E-17176 Sant Esteve de Bas. GI Spain, 
spectrophotometer: Jenway model 6405 UV/Vis). In addition, 
hemoglobin (Hb) level was measured in 0.04 mL of whole blood 
diluted in Drabkin reagent (Hycel de Mexico SA de CV), the 
hematocrit (Ht) percentage was determined by capillar)- tube 
centrifugation (10000 rpm for 5 min), and the glucose (Glc) 
concentration was measured from a 0.6 |.lL capillary blood sample 
using a portable glucometer (Bayer, CONTOUR TS, No Coding, 
Ascensia MICROLET, Japan). 

Maximal Oxygen Uptake Value 

V02max was determined using a standardized maximum 
exercise effort protocol on a treadmill, according to the 
recommendations of the American College of Sport Medicine 
[30]. A Quark b2. 8.0 bV oximeter (Cosmed, Rome, Italy) was 
used to measure oxygen (O2) consumption and carbon dioxide 
(CO2) exhalation during testing. Each participant performed a 
simulation of the test for 4-7 days prior to the stress test, so that 
they were aware of the effort required. Participants were advised 
to avoid consuming alcohol and caffeine as well as avoid smoking 
and vigorous physical activity for at least 48 h prior to testing. The 
stress test was performed 3-4 h after food intake. 

The Aerobic Load 

This consisted of a 30-mrn exercise on a treadmill, followed by 
recovery in the supine position for 30 min. To achieve greater 
precision in load development, the test was structured with five 
periods as follows: (a) Rest, 6-10 min sitting at rest; (b) Warm-up, 
5-min running on the treadmill to gradually reach 40% of 
V02max; (c) Aerobic exercise, 30 min at 60% of V02max; (d) 
Cool-down, 2 min of gradually reducing the running speed to 
zero; and (e) Post-exercise recovery, 30 min at rest in the supine 
position. (Figure 1). 

The CO2 production (VCO2) and O2 uptake (VO2) were 
measured throughout the stress test. RQnp, absolute oxygen 
consumption (VO2 in mL/min), relative oxygen consumption 
(VO2 in mL/kg/min), and relative CO2 production (VCO2 in 
mL/Kg/min) were calculated during all 5 aerobic load periods. 

Venous blood samples were collected from the cubital vein at 4 
times: sample 1, at 3—4 min before the warm-up; sample 2, at 
2 min post-exercise recovery; sample 3, at 15 min post-exercise 
recovery; and sample 4: at 30 min post-exercise recovery 
(Figure 1). Blood samples were used to determine levels of the 
following: epinephrine (Epn), norepinephrine (Nep), insulin, 
glucagon, Cortisol (Cor), growth hormone (GH), Gly, FFAs, 
beta-hydroxybutyrate (beta-HB), Glc, lactate (Lac), albumin 
(ALB), and Ht. Plasma hormone level measurements were 
performed by the UNIDAD DE PATOLOGIA CLINICA 
(Mexico Av. 2341. Guadalajara, Mexico), a laboratory certified 
by The College of American Pathologist (AV-ID: 1200227. LAP: 
5832901). Glc and Ht ("'0) k-vels were measured as described in the 
blood chemistry section. Lac levels were measured using a portable 
lactometer (Accutrend Plus System, Roche Diagnostics, San Cugat 
del Valles, Barcelona Spain). Standard colorimetric methods were 
used to determine beta-HB, Gly, and FFA levels (Randox 
Laboratories Ltd, Ardmore, Diamond Road, Crumlin, Co. 



Antrim, UK, BT29 4QY) and ALB levels (SAU Spinreact Ctra 
Santa Coloma, 7E-17176 Sant Esteve de Bas, GI Spain). 

The aerobic challenge was performed 10—20 days after 

maximum exercise testing to determine V02max. Women 
participants performed the aerobic load test within the first 10 
days of their menstrual cycle. 

Blood Hormone Levels 

Samples were taken once before the aerobic load challenge and 
thrice during recovery (2, 15, and 30 min in post-exercise recovery 
period) (Figure 1). At each sampling point, 16 mL of venous blood 
was obtained and used as follows: (a) 8 mL in 2 tubes (4 mL each) 
with sodium heparin to determine Epn and Nep by chemilumi- 
nescence assay, (b) 4 mL in cold EDTA for determining plasma 
glucagon levels by radioimmunoassay, (c) 4 mL in a dry tube for 
determining serum insulin and Cor levels by electrochemilumines- 
cence immunoassay and GH by chemUuminescence assay. 

Statistical Analyses 

We used SPSS version 15 (IBM SPSS Software, USA) for all 
analyses. The effects of E53 IK and W64R on the oximetry values 
at the end of each of the 16-measurement stages and blood 
hormone levels were tested both separately and together 
(interaction). The normality and variance homogeneity of each 
distribution were determined by the Kolmogorov-Smirnov and 
Levene tests, respectively. The values are presented in tables as 
mean ± SD. The mean values at each stage and time point were 
compared between K allele carriers vs. K allele non-carriers (K-l- 
and K— , respectively), R allele carriers vs. R allele non-carriers 
(R+ and R— , respectively), and, carriers of both (K+/R+), vs. 
carriers of K only (K+/R— ) or R only (R+/K— ). For normal and 
homogeneous distributions, one-way analysis of variance and post- 
hoc Tanhame tests for pairwise comparisons were used. Other- 
wise, the Kruskal-Wallis test was used to compare 3 distributions, 
and the Mann-Whimey U-test for pairwise comparisons. The 
statistical significance level was set at/)<0.05. 

Results 

Participant Characteristics 

Of the 45 individuals selected, only 37 (18 men, 19 women) 
successfully completed the experiment; 3 were excluded because of 
heart problems and 5 were lost to follow up (they relocated before 
the study was completed). All parents and grandparents were born 
in Mexico, with 35 originating from Western Mexico and two 
from South-Central Mexico. (Table SI in File SI). The clinical 
characteristics were as follows: age range, 22.9±3.1 years; height, 
162.7±7.8cm; weight, 63.2±10.4kg; BMI, 23.8±3.1 kg/m^; 
systolic blood pressure, 115.49±13.4 mmHg; and diastolic blood 
pressure, 71.6±6.9 mmHg. 

Data Sorting 

AU functional measurements were grouped and compared 
according to the genetic variant studied, as follows: (a) R+ (20; 10 
males, 10 females) versus (vs.) R— (17; 8 males, 9 females) to assess 
the effect of WGIR: ib) K+ (27: 13 males, 14 females) vs. K- (10; 5 
males, 5 females) to assess the c-ffect of E531K; and (c) K+/R+ (12; 
6 males, 6 females) vs. K+/R— (15; 7 males, 8 females) vs. K— / 
R+ (8; 4 males, 4 females) to test for an interaction between alleles. 
Only two participants (1 male and 1 female) were K— /R— ; 
therefore, this genotype was not included in the statistical 
comparison (Table S2 in File SI). 
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Figure 1. Aerobic load structure (5 exercise periods and 4 blood sample Times) and exercise intensity programme. 

doi:10.1371/journal.pone.0096791.g001 



Anthropometric, Dietary, and Blood Characteristics 

There were no significant differences in age or adiposity indices 
(BMl, WC, S9S, S4S, and WHR) between genotypes within the 3 
comparison groups (R-l- vs. R— ; K-l- vs. K— ; K-I-/R-I- vs. K-I-/R— 
vs. K— /R-I-; ^>0.05) (Table 1). Furthermore, the mean dietary 
protein, fat, carbohydrate, dietary SFA/kg, and MUFA/kg intakes 
did not significantly differ between groups. The amount of PUFA/ 
kg was marginally higher for K— than for K-l- (/)<0.070) (Table 2). 
Measured levels of Glc, Ht, Hb, TAG, T-Chol, HDL-C, and 
LDL-C were not significantiy different between genotypes within 
the comparison groups (/)>0.05), indicating that basal metabolic 
status was similar (Table 3). The mean HR values at 1 min after 
exercise in the Tecumseh step test did not significandy differ 



between any groups (Table S3 in File SI), establishing similar 
levels of cardiopulmonary fitness between the groups. 

Oximetry Indicators 

The VOamax percentage ("/oVOjmax) was taken as the main 
indicator of exercise intensity. Participants performed the warm-up 
phase to 35.5% ±5.3% of V02max for 5 min, the aerobic exercise 
phase at 62.3% ± 1. 1% of V02max for 30 min, and the cool-down 
phase at 42.2% ±4.6% V02max for 2 min. A significant difference 
was detected between R-l- and R— groups (63.1% vs. 61.3%; 
/) = 0.027) at 10 min into the aerobic exercise phase (Table 4). 
Despite similarities in baseline metabolic values, fitness levels, 
exercise intensity, and post-exercise hormone levels (below), there 
were significant differences in RQiip values between genotypes. 



Table 1. Mean (SD) values of age, BMl and Indicators of subcutaneous fat in each of the genetic classifications. 





Genetic Classification 


Groups 


Age (years) 


BMl (kg/m^) 


WC (cm) 


^95 (mm) 


j;4S (mm) 


CPTl B alleles 


K+ 


23.1 (3.4) 


24.2 (3.2) 


77.3 (8.1) 


150.7 (50.5) 


60.2 (22.0) 




K- 


22.5 (2.1) 


22.7 (2.7) 


75.7 (8.3) 


143.9 (50.5) 


55.5 (21.1) 






"p = 0.663 


"p = 0. 182 


''p = 0.618 


^p = 0.720 


"p = 0.562 


ADRB3 alleles 


R+ 


22.4 (1.9) 


23.6 (3.7) 


76.8 (8.5) 


148.2 (53.1) 


57.5 (22.7) 




R- 


23.5 (4.1) 


24.0 (2.4) 


76.8 (7.8) 


149.6 (47.4) 


60.6 (20.7) 






"p = 0.266 


"p = 0.659 


''p = 0.996 


'^p= 0.935 


"p = 0.668 


Interaction ADRB3-CPT1 B 


K+m+ 


22.3 (1.9) 


24.8 (4.1) 


79.0 (8.5) 


159.6 (55.8) 


61.1 (23.5) 




K+m- 


23.6 (4.3) 


23.7 (2.2) 


75.8 (7.8) 


143.5 (46.4) 


59.5 (21.5) 




R-I-/K- 


22.5 (2.1) 


21.7 (1.8) 


73.6 (7.8) 


131.1 (46.9) 


52.2 (21.9) 






''-"'p = 0.973 


''-"'p = O.OSI 


''-"'p = a;92 


''-"'p = 0.402 


''-'^p = 0.669 



SD: Standard deviation. A: One-way analysis of variance (ANOVA), significance p<0.05. K-W. Kruskal-Wallis test, significance p<0.05. 

BMl: Body mass index. WC: Waist circumference. ^9S: Sum of triceps, subscapular, biceps, iliac crest, supraspinal, abdominal, anterior thigh, medial calf, and forearm 
skin folds. ^4S: Sum of subscapular, iliac crest, supraspinal, and abdominal skin folds. {Skinfold of trunk). K+: K allele carriers. K— : K allele non-carriers. R+: R allele 
carriers. R— : R allele non-carriers. K+/R+: K and R allele carriers. K+/R— : K allele only carriers. R+/K— : R allele only carriers, 
doirl 0.1 371 /journal.pone.0096791 .tOOl 
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specifically between K+ and K— groups (Figure 2A) with either R+ 
or R— (Figure 2C). The RQnp values were not significantly difierent 
between genotypes within the comparison groups for the rest, warm- 
up, aerobic exercise, and cool-down periods, but the RQnp value in 
the K+ group remained significandy higher at 10 (/) = 0.003), 20 
(/) = 0.003), and 25 min (p = 0.048) after exercise when compared 
with the K— group. Moreover, the RQnp value was significandy 
higher in the R— than in the R-l- participants at 30 min after exercise 
(p = 0.027) (Figure 2B). The R+/K- group had an RQnp value that 
was significantly lower than that of the K+/R+ and K+/R— groups 
at 10 min after exercise (p = 0.022), and significantly lower than that 
of the K+/R— group at 20 min after exercise (p - 0.028) (Figure 2C) 
(Table 5). Mean HR were no significandy different between 
genotypes throughout the exercise and recovery periods {p>0.05) 
(data not shown). 

Blood Indicators of Aerobic Load 

Lac levels were not significantly difierent between genotypes 
within the 3 comparison groups prior to testing and at 2, 15, and 
30 min after exercise (/)>0.05). Lac levels exceeded 2 mM only at 
2 min after exercise, and lower values were measured at 15 and 
30 min. Glc levels decreased slightly during the post-exercise 
r(;c()\'er\' phase, but not significant difierences were observed 
between the genotypes. Similarly, ALB levels and Ht were not 
significandy difierent between genotypes at any stage (Table S4 in 
File SI). The K— group had significantly higher FFA levels than 
those in the K+ group at 2 min after exercise (J) = 0.044) (Figure 3B; 
Table 6), whereas there were no significant difiFerences in Gly levels 
between genotypes at any of the 4 measurement points. Beta-HB 
levels tended to increase during the post-exercise recovery phase 
(Table 6); however, no significant differences were found between 
the R+ and R— participants (Figure 4A), and levels were 
significandy higher at 15 min (/; = 0.022) and 30 min (/> = 0.008) 
after exercise in the K— than in the K+ participants (Figure 4B). 
Furthermore, a significantly higher beta-HB level was detected at 
30 min after exercise in the K— /R+ group compared with either the 
K-I-/R-I- or K-I-/R- groups (p = 0.029) (Figure 4C). 

Plasma Hormone Levels during Aerobic Load 

Insulin, Cor, and GH levels did not differ significandy between 
genotypes within the comparison groups at any time. Glucagon 
levels remained below 50 pg/ mL at all 4 measurement times. Epn 
and Nep levels reached peak levels at 2 min after exercise (84.81 
and 676. 14 pg/ mL, respectively) and declined thereafter. The only 
significant difference in Nep level was detected at 30 min after 
exercise between R+ and R— participants (230.6 pg/mL vs. 301.7 
pg/mL; p = 0.028) (Table S5 in FUe SI). In general, the hormonal 
responses were as expected and were generally similar among the 
genotypes. 

Discussion 

The increasing rate of fatty acid oxidation induced by aerobic 
exercise is modulated by multiple factors. In this study, we 
therefore sought to reduce the influence of other polymorphic 
genes involved in the lipolysis and fatty acid oxidation pathway by 
including only Mexicans from only 2 generations, and primarily 
from the same region of Western Mexico. 

Carriers of the K allele exhibited a higher mean RQnp during 
the post-exercise recovery phase, indicating reduced lipid beta- 
oxidation compared with non-carriers. In light of previous studies 
demonstrating an association between the K allele and obesity 
indices [1 1], this study suggests that the K allele may predispose to 
decreased FFA oxidation when the energy demand from lipids is 



increased during the recovery from effort. This can be a factor 
contributing to the accumulation of fat in adipose tissue given 
similar dietary and activity conditions to non-carriers. 

The difference in RQnp between K allele carriers and non- 
carriers cannot be explained by differences in diet, general 
physical condition, or demographic valuables (sex ratio or age). 
Expression of the human CPTIB gene in transgenic mice 
increased 2-3-fold during both fasting and high-fat diets [31], 
whereas PUFA modulate the expression of genes involved in lipid 
metabolism, such as the peroxisome proliferator artivated receptor 
[32] . Thus, dietary factors can inffuence the expression of CPTIB 
and other genes involved in fat mobilization. However, the 
comparison groups reported similar dietary macronutrient com- 
positions (including SFA, MUFA, and PUFA). In addition, the 
level of cardiopulmonary fitness did not differ between genotypes 
(Table S3 in File SI); thus, obviating the differential effects of 
adaptation to physical activity [33,34]. Moreover, subcutaneous 
fat, Glc, Hb, lipid, hpoprotein, and Ht levels did not differ between 
the comparison groups, indicating that metabolic conditions were 
similar at basehne. During the development of aerobic load, there 
were no differences in %V02max among the genotypes. 
Moreover, blood Lac levels did not reach 4 mM at any point 
after exercise, which is indicative of aerobic metabolism (Table S4 
in File SI). These results confirmed that exercise load was 
predominantiy aerobic and of similar intensity among all 
participants. 

Hormone levels involved in the modulation of lipolysis and fat 
oxidation did not differ significandy among the genotypes after 
exercise. Furthermore, ALB Hb, and Ht levels were not 
significandy diflFerent. The apparent differences in lipid metabohsm 
between K allele carriers and non-carriers was therefore unlikely to 
be due to differences in hormonal patterns, basal FFA, or oxygen 
transport and FFA transport capacity to active muscle tissue. 

No significant differences were observed in plasma FFA and Gly 
levels between R allele carriers and non-carriers. Moreover, the 
similarity in catecholamine and other hormone levels before and 
after exercise indicates that the R allele does not affect the 
activation of lipolysis; these genot)ipes probably had similar FFA 
metabolism rates. There was a trend toward high values among 
R+ in RQnp, but only reached significant difference at 30 minute 
post-exercise recovery period compared to non-carriers (possibly 
related to the 12 K allele carriers) and a downward trend among 
non-carriers (possibly related to the 15 E allele carriers). In 
addition, it has been observed that LCFA oxidation in the 
mitochondria increases with FFA availability [9,35], which was 
similar in both groups. In addition, beta-HB levels did not differ 
betwcx-n R allele carriers and non-carriers during recovery. These 
results indicate that the availability of FFA, mobilized from 
visceral adipose tissue was similar to its consumption by skeletal 
muscle tissue. 

Most participants in this study were heterozygous for the R 
allele, which already been reported in other studies [10,14,16,18]. 
Thus, the gene dose may have contributed to the similarity in 
Kpolysis indicators between R allele carriers and non-carriers. In 
addition, these indices were measured for only 30 min after 
exercise; delayed responses may therefore exist with measureable 
differences. Another possibility is that the effect of the W64R 
variant of ADRB3 could be masked by the effect of other lipolysis 
agents, including cytokines that are activated during muscle 
contraction [36]. Nonetheless, several studies and meta-analyses 
have shown that although the R allele is associated with increased 
BMI and insulin resistance, its contribution to the phenotype is 
weak [14,17,38]. Our results indicate that the W64R allelic variant 
ot ADRB3 does not affect the activation of lipolysis. 
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Figure 2. Non-protein respiratory quotient (RQnp) values during aerobic load. (A) R allele carriers (1 2 K and 8 E carriers) versus R allele non- 
carriers (15 K and 2 E carriers). (B) K allele carriers (12 R and 15 W allele carriers) versus K allele non-carriers (2 R and 8 W carriers) The asterisk and 
annpersand indicate a significant difference. (C) Interaction (CPT1B-ADRB3). The asterisk and ampersand indicate the significant difference between 
the R allele carrier group, with 8 carriers of the E allele, and the other groups. 
doi:1 0.1 371/journal.pone.0096791 .g002 
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Figure 3. Free fatty acid (FFA) levels before and after tKie aerobic load. (A) R allele carriers (12 K and 8 E carriers) versus R allele non-carriers 
(15 K and 2 E carriers). (B) K allele carriers (12 R and 15 W allele carriers) versus K allele non-carriers (2 R and 8 W carriers). The asterisk indicates 
significant difference. (C) Interaction {CPTl B-ADRB3), carriers of both K and R alleles versus only K and only R allele carriers. Venous blood was 
sampled at rest and at 2, 15, and 30 min into the recovery period. 
doi:l 0.1 371/journal.pone.0096791 .gOOB 
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Table 6. Mean free fatty acid and beta-hydroxy butyrate at rest and at 2, 15, and 30 min post-exercise for each of the genetic 
classifications. 





Genetic Classification 


Rest 


2' post 


15' post 


30' post 






exercise 


exercise 


exercise 


FFA (ntiM) K+ 


0.24 (0.21) 


0.82 (0.51)* 


0.82 (0.51) 


0.65 (0.45) 


K- 


0.32 (0.19) 


1.16 (0.53)* 


1 .09 (0.46) 


96.2 (0.38) 




"p = 0.067 


"p = 0.044 


"p = 0.!J9 


"p = 0.062 


R+ 


0.26 (0.17) 


0.91 (0.48) 


0.88 (0.42) 


0.74 (0.39) 


R- 


0.26 (0.32) 


0.92 (0.59) 


0.90 (0.60) 


0.70 (0.51) 




'*p= 0.972 


''p = 0.954 


'^p = 0.901 


"p = 0.785 


K+/R+ 


0.20 (0.13) 


0.69 (0.34) 


0.74 (0.35) 


0.60 (0.33) 


K+/R- 


0.27 (0.35) 


0.91 (0.61) 


0.88 (0.61) 


0.69 (0.53) 


R+/K- 


0.35 (0.19) 


1.22 (0.51) 


1.10 (0.45) 


0.95 (0.40) 




''p= 0.627 


"^p = 0.189 


"p = 0.464 


''p = 0.414 


beta-HB (mM) K+ 


0.079 (0.05) 


0.089 (0.06) 


0.167 (0.09)* 


0.189 (0.11)* 


K- 


0.070 (0.05) 


0.090 (0.06) 


0.240 (0.08)* 


0.300 (0.09)* 




"p = 0.699 


"p = 0.950 


"p = 0.022 


"p = 0.008 


R+ 


0.070 (0.05) 


0.085 (0.06) 


0.185 (0.09) 


0.220 (0.12) 


R- 


0.082 (0.05) 


0.094 (0.06) 


0.188 (0.10) 


0.218 (0.11) 




"p = 0.491 


"p = 0.618 


"p = 0.961 


"p = 0.987 


K+/R+ 


0.075 (0.05) 


0.075 (0.06) 


0.150 (0.08) 


0.167 (0.11) 


K+/R- 


0.080 (0.06) 


0.100 (0.05) 


0.180 (0.10) 


0.207 (0.11) 


R+/K- 


0.062 (0.05) 


0.100 (0.05) 


0.240 (0.09) 


0.300 (0.09)** 




"p = 0.756 


''p = 0.444 


"p = 0.098 


''p = 0.029 



SD: Standard deviation. A One-way analysis of variance (ANOVA). U Mann-Whitney test. K: Kruskal-Wallis test. Significance p<0.05. Asterisk (•) indicates the values that 
significantly differ in the same stage. Double asterisk (**) indicates that the value differs significantly from the other values in the same stage. 
FFA: Free fatty acids. beta-HB: Beta hydroxybutyrate. K+: K allele carriers. K— : K allele non-carriers. R+: R allele carriers. R— : R allele non-carriers. 
doi:l 0.1 371 /journal.pone.0096791 .t006 



The RQnp values during the post-exercise recovery phase were 
higher in K allele carriers (Figure 1 -B) and independent of the R 
allele (Figure 1-C). This is consistent with the hypothesis that the K 
allele affects the translocation capacity of CPTIB, and that the effect 
is manifested under increased energy demands from fat that occur at 
the end of a prolonged workout. K allele carriers showed higher 
RQnp values during recovery, suggesting contributions to energy 
generation from fat equivalents of 19.3% at 10 min, 49.3% at 
15 min, and 63.1% at 25 min compared witli 59.7%, 88.0%, and 
91.6% at the same post-exercise times in non-carriers [4]. Thus, K 
allele carriers demonstrate significantly lower fat metabolism than 
non-carriers. It is unlikely that the higher RQnp values among K 
allele carriers resulted from a reduced availability of circulating FFA 
because there were no differences in plasma FFA concentrations 
between groups after exercise, in addition the Gly values were 
similar, and there were no differences in hormone levels or exercise 
intensity. However, this is an opportunity for future research. Napal 
et al. [39] determined that the in vitro substitution of glutamic acid 
at position 590 of CPT 1 A increased sensitivity to malonyl coenzyme 
A. They concluded that the glutamic acid residue may impart 
structural stability to the protein, and thus affect the functional 
properties of the whole CPT family. In this study, we found that the 
substitution of glutamic acid at position 53 1 by lysine in the CPTIB 
protein had a significant effect on RQnp increase. 

Notably, the sample size of this study was too small to include 
sufficient numbers with the rare K— /R— genotype; therefore, we 
could only assess the effects of allelic interaction among carriers of 
both alleles (K-I-/R-I-) and those carrying only one of the alleles 



(K-H/R- and K-/R-H). Consequently, the K-/R-H group also 
functioned as the non-carrier group for the K allele. During the 
recovery period, FFA levels were higher in the K— /R-l- group 
than in the other two groups, although the difference was not 
statistically significant. At 1 0 and 20 min after exercise, this group 
also showed significandy lower RQnp values (Figure 2-C), and at 
30 min after exercise showed significantly higher beta-HB levels 
(Figure 4-C). These observations indicate that the E531K had 
significant effects that were independent of W64R. 

In this experimental design, the observation of the influence of 
the E531K SNP of CPTIB was facihtated because the CPTl 
protein represents the only route of entry of LCFA into the beta- 
oxidation pathway. In contrast, ADRB3 lipolysis is modulated by 
the receptor as well as by several hormonal and metabolic factors. 
In this study, we ensured only that there were no significant 
differences in plasma levels of catecholamines, insulin, glucagon. 
Cor and growth hormone, lipids, diet (particularly PUFA content), 
exercise intensity, body adiposity, and cardiorespiratory fitness. 
Therefore, further studies are needed to strengthen and clarify the 
functional impact of the gene mutations involved in lipid 
metabolism, together with other factors such as the levels of 
interleukiii-6 produced by muscle during contraction that can 
induce lipolysis within adipose tissue that peaks at the end of 
exercise [32,37] To this end, we have estabhshed a DNA bank 
with 285 participants to facilitate future cohort studies. 
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Figure 4. beta-Hydroxybutyrate (beta-HB) levels before and following the aerobic load. (A) R allele carriers (1 2 K and 8 E carriers) versus R 
allele non-carriers (15 K and 2 E carriers). (B) K allele carriers (12 R and 15 W allele carriers) versus K allele non-carriers (2 R and 8 W carriers). The 
asterisk indicates significant difference. (C) Interaction (CPT1B-ADRB3), carriers of both K and R alleles versus only K and only R allele carriers. The 
asterisk indicates a significant difference between the R allele group with 8 carriers of the E allele and the other groups. Venous blood was sampled at 
rest and at 2, 15, and 30 min into the recovery period. 
doi:l 0.1 371/journal.pone.0096791 .g004 
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Conclusion 

The substitution of glutamic acid at position 531 by lysine in the 
CPTIB protein decreases the mitochondrial beta-oxidation 
pathway, which increases the RQnp value during recovery from 
exercise. This may in turn contribute to either weight gain or a 
reduced weight-loss following exercise. 

Supporting information 

FUe SI Contains Tables S1-S5. Table SI: Birthplace of the 
participants who completed the experimental tests. Table S2: 
Participants distribution in each comparison groups. Table S3: 
Mean (SD) values of heart rate (HR) in Tecumseh Step Test. 
Table S4: Means values and (SD) of blood indicators according 
with genetic classification. Table S5: Means and (SD) of blood 
hormones concentrations and hematocrit according with genetic 
classification. 
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